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Abstract.  The  ion  velocity  distribution  (IVD)  is  important  in  plasma  etching  of  microfeatures.  IVD  at  a  rf 
biased  wafer  is  studied,  first  analytically  using  probability  theory  and  then  numerically  by  using  a  particle 
simulation  method.  The  analytic  expression  shows  that  IVD  is  governed  by  the  parameter  qVrt where 
q  is  the  charge  of  ion,  Kf  is  the  rf  bias  amplitude,  cu  is  the  rf  bias  angular  frequency,  I  is  the  penetration 
depth  of  bias  potential,  and  m;  is  the  mass  of  ion.  The  analytical  expression  is  applicable  to  the  case  when 
the  ion  collisions  in  the  penetration  depth  are  negligibly  few  and  the  rf  period  of  biasing  is  much  shorter 
than  the  time  that  ions  take  in  traversing  the  depth  1.  The  IVDs  for  general  conditions  are  also  examined 
using  the  self-consistent  particle-in-cell/Monte  Carlo  simulation. 


I  INTRODUCTION 

Weakly  ionized  plasmas  used  in  materials  processing  generally  consist  of  positive  and  negative  ions,  electrons, 
and  neutral  species.  In  the  case  when  electrons  are  more  abundant  than  negative  ions,  a  negative  potential 
is  formed  near  the  surface  of  the  floating  wafer  immersed  in  the  plasma.  In  processing  plasmas  the  electron 
temperature  is  few  electron  volts,  while  the  ion  temperature  is  roughly  equal  to  the  background  gas  temperature. 
To  keep  the  balance  between  the  negative  and  positive  currents  on  the  floating  surface,  low  energy  electrons 
should  be  reflected  in  the  sheath  region  back  to  the  plasma  bulk.  This  is  the  physical  reason  why  the  potential 
of  the  floating  wafer  is  lower  than  the  plasma  potential. 

Negative  sheath  potential  (measured  from  the  plasma  bulk)  accelerates  the  positive  ions  to  the  wafer.  If 
no  collision  is  assumed  in  the  sheath  region,  the  positive  ions  impinge  almost  perpendicularly  onto  the  wafer. 
However,  electrons  are  decelerated  and  lose  the  velocity  component  normal  to  the  wafer. 

When  ions  and  electrons  come  into  a  microscopic  trench,  electrons  are  easily  sticked  near  the  top  of  the 
side  of  the  trench  wall  while  ions  can  reach  the  bottom  of  the  trench.  This  kind  of  charge  separation  inside 
microfeatures  causes  a  charge  build-up.  The  strongly  localized  electric  field  due  to  the  charge  build-up  distorts 
the  trajectories  of  incident  ions,  which  results  in  an  anomalous  etch  profile  called  “notching”.  The  localized 
charge  build-up  also  causes  the  electrical  breakdown  induced  by  fatal  tunneling  currents  through  gate-oxides. 

In  order  to  suppress  the  so-called  charging  damages,  new  ideas  on  plasma  sources  are  proposed,  e.g.,  time- 
modulation  of  the  source  power  [1]  and  use  of  ultra-high  frequency  (500  MHz)  power  [2] .  These  approaches  are 
intended  to  lower  the  electron  temperature  and  hence  generate  more  negative  ions.  To  understand  the  role  of 
negative  ions,  let  us  consider  the  plasma  including  only  positive  and  negative  ions;  when  the  wafer  is  biased 
at  low  frequency  (<  3 MHz),  the  positive  and  negative  ions  impinge  on  the  bottom  of  the  trench  alternately. 
Thus  no  charge  build-up  occurs. 

The  ion  energy  distribution  (lED)  and  ion  angle  distribution  (lAD)  of  ions  that  impinge  on  a  biased  wafer 
are  important  to  clarify  the  mechanism  of  etching  and  charging  damage  of  microfeatures.  In  low  gas  pressure 
and  high  plasma  density,  the  rf  wafer  biasing  is  employed  for  controlling  incident  ion  energy  independently  of 
the  source  power.  Sobolewski  et  al.  [3]  measured  the  lEDs  for  various  bias  amplitudes  and  frequencies  in  a 
planar  inductively-coupled  plasma  (ICP)  reactor.  Hoekstra  and  Kushner  [4]  calculated  the  lEDs  in  an  ICP. 
They  discussed  the  effects  of  the  bias  amplitude,  gas  pressure,  and  source  power  on  the  lED.  Notching  can  be 
explained  from  charging  of  the  surfaces  of  microfeatures.  Kinoshita  et  al.  [5]  presented  the  relationship  between 
the  charging  potential  and  the  lED.  Hasegawa  et  al.  [6]  reported  that  an  appropriate  rf  biasing  suppresses  the 
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charging  damage. 

Nasser  et  al.  [7]  found  that  the  penetration  depth  of  the  potential  oscillation  due  to  rf  biasing  is  much  larger 
than  the  Debye  length.  Thus  the  lED  is  greatly  affected  by  the  potential  oscillation.  In  this  report  we  consider 
two  approaches  for  obtaining  the  IVD  at  a  rf  biased  wafer.  First,  using  the  probability  theory,  we  obtain  the 
IVD  analytically  under  certain  simplified  situations.  The  obtained  expression  shows  that  the  IVD  is  governed 
by  a  similarity  plasma  parameter.  Second,  using  self-consistent  particle- in-cell/Monte  Carlo  simulation,  we 
determine  the  oscillatory  sheath  potential  due  to  rf  biasing  and  hence  obtain  the  IVDs. 


II  ANALYTICAL  ION  VELOCITY  DISTRIBUTION 


First  we  obtain  the  IVD  analytically  using  the  probability  theory  [8].  Ahn  [9]  measured  the  instantaneous 
plasma  potential  in  the  ICP.  The  potential  that  is  time-modulated  by  high  frequency  biasing  is  well  approxi¬ 
mated  by  the  exponentially  attenuating  sine  curve.  In  this  case,  the  velocity  of  ion  normal  to  the  wafer, 
can  be  expressed  as  =  Uw  -fUw,  where  is  the  contribution  from  a  time-averaged  potential  and  is  that 
from  an  oscillatory  part  of  the  potential  [10].  For  high  frequency  biasing,  we  found  that  the  time  when  the 
ion  impinges  on  the  wafer  is  uniformly  distributed.  Assuming  no  collision,  we  can  then  obtain  the  probability 
density  function  G{v^)  for  the  time-modulated  velocity 


G{v^) 
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(1) 


where  A  is  given  by 
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qVti 
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Here  mi  is  the  mass  of  ion,  q  is  the  charge  of  ion,  Uf  is  the  amplitude  of  the  rf  bias,  I  is  the  penetration  depth 
of  the  biasing  potential,  and  ui  is  the  bias  angular  frequency.  The  probability  density  function  for  tJw 

can  be  obtained  by  assuming  that  ions  in  the  plasma  bulk  are  in  equilibrium.  Now  the  probability  density 
function  F||(uw)  for  the  sum  Uw  =  Uw  +  Vy,  is  given  by  the  convolution  of  G(uw)  and  H{vy,) 

■P’||(^w)  =  j  G{v^ -v)H{v)dv.  (3) 

To  check  the  validity  of  the  analytical  IVD,  we  separately  determined  the  IVD  from  numerical  solutions  of  the 
equation  of  motion  and  compared  the  two  results.  The  best  fitting  of  Ahn’s  [9]  experimental  data  of  potential 
becomes  t)  =  4.5  +  0.019a;  -f  2.4  x  10“^x^  —  16exp(— x/2.1)  +  55exp(— x/2.8)  sin(a;t),  where  x  in  units  of 
mm  is  the  normal  distance  measured  from  the  wafer  at  x  =  0.  The  rf  bias  frequency  is  15  MHz.  We  calculated 
the  motion  of  ClJ  ions.  Figure  1  shows  the  comparison  of  the  analytical  and  numerical  IVDs.  The  analytical 
expression  shows  a  good  agreement  with  the  numerical  solution.  In  Fig.  I,  the  distance  between  the  two  peaks 
is  2A(=  2qVri/miUjl). 


Ill  SELF-CONSISTENT  PIC/MC  SIMULATION 

It  is  clear  that  the  IVD  is  a  functional  of  the  sheath  potential  ^(x,  t),  which  depends  on  the  wafer  biasing  and 
collisions.  In  the  case  when  there  is  no  experimental  data  on  (j){x,  t),  it  is  necessary  to  calculate  (f){x,  t)  and  IVD 
simultaneously.  In  order  to  examine  the  effects  of  bias  frequency  and  ion  collisions  on  IVD,  the  self-consistent 
particle- in-cell/Monte  Carlo  simulation  is  employed. 

The  background  gas  is  argon  which  is  assumed  to  be  uniform  and  in  equilibrium  at  temperature  of  300  K. 
Although  e“,  Ar+,  and  Ar  are  taken  into  consideration  in  the  present  calculation,  we  have  only  to  examine 
the  motion  of  e“  and  Ar"*" .  The  velocity  of  Ar  is  necessary  only  when  determining  the  postcollisional  velocity 
of  e-  or  Ar+. 

Since  the  number  density  of  charged  species  is  less  than  10^'^  m“^  in  the  present  problem,  Coulomb  collisions 
(e“-e“,  Ar^-Ar"*"  and  e^-Ar"*")  can  be  neglected  and  therefore  we  consider  e“-Ar  and  Ar^-Ar  collisions.  For 
e“-Ar  collisions,  elastic  collision,  ionizing  collision,  and  also  25  exciting  collisions  are  taken  into  consideration 
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FIGURE  1.  Ion  velocity  distributions  (no  collision). 


Electron  Energy  e  (eV) 
FIGURE  2.  Cross  sections  for  e“-Ar  collision. 


based  on  a  set  of  the  cross  section  data  of  Kosaki  and  Hayashi  [11].  As  for  ionization,  Peterson  and  Allen’s  [12] 
integral  cross  section  is  multiplied  by  0.9  to  match  it  with  Kosaki  and  Hayashi’s  data.  The  total  number  of 
collisional  events  is  27.  Each  event  is  expressed  hy  k  {k  =  1,2,  •••,27);  k  =  1-25  for  excitation,  fc  =  26  for 
ionization,  and  k  =  27  for  elastic  collision.  The  cross  section  for  e“-Ar  collision  is  shown  in  Fig.  2.  The 
probability  of  occurrence  of  the  kth  event  in  time  At  is 

P,{k)  =  NB<7k{e)(^)  At,  (4) 

where  Nb  is  the  number  density  of  background  argon  gas,  e  is  the  electron  energy,  and  me  is  the  electron  mass. 

Elastic  collision  and  resonant  charge  exchange  are  considered  for  the  Ar+-Ar  collision.  The  simple  and 
computationally  efficient  model  proposed  in  our  previous  work  [13],  the  use  of  which  was  shown  to  reproduce 
the  measured  data  of  ion  drift  velocity  and  transverse  diffusion  coefficient,  is  employed.  The  model  is  a 
combination  of  isotropic  scattering  and  charge  switching,  both  occurring  with  the  same  probability.  The 
probability  Pi  that  an  ion  Ar+  collides  with  atom  Ar  in  time  At  is  given  by 

=  NB9Txd\BAt,  (5) 

where  g  is  the  relative  velocity  between  Ar+  and  Ar  averaged  over  the  Maxwellian  distribution  for  Ar  with 
temperature  Tb(=  300  K)  and  is  the  collision  cross  section  twice  as  large  as  Trdig,  cIbb  being  the  viscosity 
diameter  of  Ar. 

One  time  step  of  the  simulation  consists  of  solving  the  equation  of  motion  for  e^  and  Ar+,  obtaining  the 
charge  density,  solving  the  Poisson  equation,  and  calculating  collisions.  The  time-evolution  is  calculated  until 
a  periodic  steady  state  is  obtained. 

The  one-dimensional  computational  domain  from  the  biased  wafer  {x  =  0)  to  the  plasma  bulk  (x  =  L)  is 
shown  in  Fig.  3.  The  unknowns  in  the  present  problem  are  the  dc  bias  Vuc  at  the  wafer  and  the  flow  velocity 
of  ions  at  X  =  L.  The  electron  flow  velocity  at  x  =  A  is  equal  to  the  ion  flow  velocity  because  of  no  net  current 
in  the  plasma  bulk.  The  bias  Vhc  is  determined  from  the  floating  condition;  the  net  current  during  one  period 
of  rf  biasing  is  zero  in  periodic  steady  state.  Ar+  and  e“  come  from  x  >  L  and  cross  the  presheath  edge  at 
X  =  L.  The  flow  velocity  is  determined  in  such  a  way  that  the  plasma  density  is  continuous  at  x  =  L. 

The  length  L  of  the  computational  domain  is  100  mm,  as  is  shown  in  Fig.  3.  Let  the  x-axis  be  in  the  direction 
normal  to  the  wafer  with  the  origin  at  the  wafer.  We  used  equally  spaced  grids  of  Ax  =  0.1  mm.  That  is,  the 
domain  was  divided  into  1000  cells.  Note  that  the  Debye  length  is  0.11mm  for  electron  density  of  10^®  m^^ 
and  electron  temperature  of  2eV.  The  electron  time  step  Ate  was  determined  so  as  to  satisfy  the  condition 
that  the  total  collision  probability  for  electron  is  less  than  0.03.  This  condition  also  meets  the  requirement  for 
Nanbu’s  method  [14]  that  the  maximal  probability  of  collisional  events  must  be  less  than  1/27,  where  27  is  the 
total  number  of  collisional  events.  From  the  requirement  we  chose  Afg  =  9.26  x  10^^^  s  at  the  background  gas 
pressure  pb  =  5mTorr,  Atg  =  5.29  x  10^^^  s  at  pb  =  lOmTorr,  and  Atg  =  2.65  x  10^^^  s  at  pb  =  20mTorr. 
The  ion  time  step  At,  was  chosen  to  be  At,  =  30Ate,  for  which  the  Ar+-Ar  collision  probability  is  less  than 
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FIGURE  3.  One-dimensional  computational  domain.  Distance  x  (mm) 

FIGURE  4.  Electric  potential  (j>.  Wafer  is  at  a:  =  0. 


0  50  100  0  50  100 


Distance  x  (mm)  Distance  x  (mm) 

FIGURE  5.  e“-Ar  elastic  collision  rate.  FIGURE  6.  Ion  flux. 

0.05.  Since  the  plasma  bulk  is  charge-neutral,  a  common  value  of  10^®  m”^  was  chosen  as  the  ion  and  electron 
densities  in  the  plasma  bulk.  In  processing  plasmas  the  electron  temperature  Tg  is  much  higher  than  the  ion 
temperature  Ti  due  to  the  power  deposition  to  electrons.  Therefore  we  set  T;  =  500  K  and  feuTe  =  2  eV  in  the 
plasma  bulk  {x  >  L),  where  fee  is  the  Boltzmann  constant. 

IV  RESULTS  AND  DISCUSSION 

The  validity  of  the  method  described  above  was  first  examined  by  calculating  the  structure  of  the  potential 
in  the  case  of  no  rf  biasing,  the  wafer  being  floating.  Figure  4  shows  the  electric  potential  (p.  The  region 
X  <  L{=  100  mm)  consists  of  the  sheath  plus  presheath.  In  the  presheath  the  charge  neutrality  is  kept.  Note 
that  the  potential  in  the  plasma  bulk  is  set  zero.  Assuming  no  collision  for  ions  in  the  sheath  and  presheath 
results  in  the  floating  potential  (pf  on  the  wafer  as  [15] 


As  is  shown  later,  the  electron  temperature  in  x  <  L  is  almost  constant  and  close  to  the  value  in  the  plasma  bulk 
(2  eV).  The  floating  potential  obtained  from  eq.  (6)  is  —10.4  V.  The  floating  potential  in  Fig.  4  is  —12.29  V  at 
the  background  gas  pressure  pb  =  5mTorr,  —12.58  V  at  ps  =  lOmTorr,  and  —12.77  V  at  pb  =  20mTorr.  The 
corresponding  flow  velocity  at  x  =  i  is  327  m/s  at  5mTorr,  112  m/s  at  lOmTorr,  and  —124  m/s  at  20mTorr. 
The  negative  flow  velocity  at  20mTorr  is  caused  by  a  considerable  increase  of  the  number  of  electrons  and 


257 


Distance  x  (mm) 
FIGURE  7.  Ion  density. 


Distance  x  (mm) 

FIGURE  8.  Number  densities  of  ion  and  electron 
(10  mTorr). 


Ion  Velocity  (km/s) 


FIGURE  9.  Ion  velocity  distribution. 
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FIGURE  10.  Electron  temperature. 


ions  in  X  <  L  due  to  ionization.  The  difference  of  (pf  between  our  data  and  eq.  (6)  is  ascribable  to  e“-Ar  and 
Ar+-Ar  collisions  in  a;  <  A.  Figure  5  shows  the  rate  for  e~-Ar  elastic  collision.  The  rate  is  proportional  to  the 
gas  pressure  and  the  electron  density  rie  approximately.  The  rate  decreases  with  x  because  tie  decreases  with 
X,  as  is  shown  later.  The  ionization  rate  is  about  4  orders  smaller  than  the  elastic  collision  rate.  As  we  have 
seen,  however,  ionization  has  a  large  effect  near  the  outer  boundary  of  the  computational  domain.  Also  ion 
(electron)  flux  is  not  constant  due  to  ionization.  This  is  shown  in  Fig.  6.  Electrons  are  decelerated  by  e“-Ar 
collisions  and  ions  are  by  Ar+-Ar  collisions.  The  fact  that  (p{  determined  from  the  simulation  is  lower  than 
the  collisionless  (pf  given  by  eq.  (6)  means  that  the  collisional  deceleration  is  stronger  for  Ar+-Ar  collision  than 
e^-Ar  collision,  thus  requiring  a  lower  (f>{  to  accelerate  ions  and  satisfy  the  floating  condition  on  the  wafer. 
The  larger  collisional  deceleration  for  Ar+-Ar  is  due  to  charge-exchange  collisions. 

Where  is  the  boundary  between  the  sheath  and  presheath?  Figure  7  shows  the  ion  density  rii.  The  electron 
density  almost  agrees  with  m  except  the  region  close  to  the  wafer.  This  region  is  enlarged  in  Fig.  8  at 
Pb  =  10  mTorr.  We  see  that  the  charge  neutrality  breaks  near  x  =  2  mm.  This  location  is  the  boundary 
between  the  sheath  and  presheath.  Figure  9  shows  the  ion  velocity  distribution  function.  The  distribution  is 
more  localized  as  pressure  and  hence,  the  ion  collision  frequency  decreases.  The  upper  limit  in  the  right  side 
of  the  peak  represents  the  ion  velocity  in  the  case  of  no  collision.  We  show  the  electron  temperature  in  Fig.  10. 
The  change  in  the  temperature  is  very  small  in  the  presheath.  As  pressure  increases,  the  electron  temperature 
in  the  presheath  slightly  decreases  due  to  the  increase  of  the  rate  of  inelastic  e^-Ar  collisions.  The  rapid 
decrease  near  the  wafer  is  due  to  the  sharp  decrease  of  the  potential  in  the  sheath.  The  ion  temperature  is  also 
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(a)  15  MHz  (b)  IMHz 

FIGURE  11.  Time-modulation  of  potential  4>  (lOmTorr). 


FIGURE  12.  Time-averaged  potential  </>  (lOmTorr). 


Distance  x  (mm) 

FIGURE  13.  Ion  density  (lOmTorr). 


almost  constant  in  the  presheath.  These  results  show  that  the  present  method  is  applicable  to  the  numerical 
analysis  of  the  sheath  and  presheath. 

We  then  applied  our  method  to  the  case  when  the  floating  wafer  is  under  rf  biasing.  We  fixed  the  gas 
pressure  at  lOmTorr  and  examined  the  effect  of  the  bias  frequency  on  the  sheath  structure.  Our  choice  is 
15  MHz  and  IMHz.  The  peak-to-peak  voltage  applied  to  the  wafer  is  fixed  at  100  V.  Figures  11(a)  and  (b) 
show  the  time-modulation  of  the  potential  (j>{x,a),  where  a{=  uit)  is  the  phase  of  the  bias  potential.  As  is 
shown  later,  the  sheath  and  presheath  boundary  is  near  x  =  5  mm.  Figure  11  shows  that  the  sheath  is  time- 
modulated  at  both  low  and  high  frequency.  Figure  11(a)  shows  that  in  case  of  high  frequency  not  only  the 
sheath  but  also  the  presheath  is  somewhat  time-modulated.  Figure  12  shows  the  time-averaged  potential  ^{x). 
The  effect  of  the  bias  frequency  on  (j)  is  small.  We  see  that  the  dc  bias  is  —59.27  V  for  15  MHz  and  —57.86  V 
for  1  MHz.  Figure  13  shows  the  ion  density.  The  higher  ion  density  at  15  MHz  is  probably  due  to  a  larger  rate 
of  ionization.  The  ion  density  is  hardly  time-modulated,  as  is  shown  in  Figs.  14(a)  and  (b).  Strictly  speaking, 
the  ion  density  is  slightly  time-modulated  at  the  low  frequency,  as  is  seen  from  Fig.  14(b).  The  electron  density 
rie  agrees  with  the  ion  density  ni  in  the  presheath.  However,  is  time-modulated  in  the  sheath,  as  is  shown 
in  Figs.  15(a)  and  (b).  The  electron  density  is  non-zero  at  the  phase  of  tt/2,  at  which  the  wafer  works  as  an 
instantaneous  anode.  The  time-averaged  electron  and  ion  densities  are  shown  in  Figs.  16(a)  and  (b).  The 
charge  neutrality  breaks  at  x  =  5  mm,  which  is  the  edge  of  the  sheath.  Although  the  plasma  density  is  higher 
at  15  MHz,  we  can  see  no  effect  of  the  bias  frequency  on  the  sheath  thickness. 
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Lastly,  the  ion  velocity  distribution  function  is  shown  in  Fig.  17.  The  gas  pressure  is  lOmTorr  as  before. 
Owing  to  Ar+-Ar  collisions,  the  function  does  not  take  the  form  of  twin  peaks  in  Fig.  1.  Since  the  highest 
velocity  is  realized  in  the  case  of  no  collision,  however,  the  right  side  of  the  peak  can  be  partly  interpreted  from 
the  analytic  expression  of  eq.  (3).  In  Fig.  1,  the  location  of  the  right  peak  is  near  to  +  A.  Since  A  in  eq.  (2) 
decreases  with  increasing  ui,  the  peak  in  Fig.  17  is  shifted  with  increasing  the  bias  frequency. 


V  CONCLUSION 

The  IVD  was  obtained  using  the  self-consistent  particle- in-cell/Monte  Carlo  simulation.  In  the  case  of  no  rf 
biasing  the  dependence  of  IVD  on  the  background  gas  pressure  was  investigated.  We  found  the  IVD  is  more 
localized  toward  the  upper  side  of  the  distribution  as  pressure  decreases.  We  then  calculated  the  oscillatory 
sheath  potential  due  to  rf  biasing.  The  IVD  under  wafer  biasing  is  broad  due  to  ion  collisions.  The  shape  of 
IVD  in  high  velocity  regime  is  analogous  to  the  analytic  IVD  for  no  collision;  the  location  of  the  peak  is  shifted 
to  a  lower  velocity  as  the  bias  frequency  increases. 
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